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We report on measurements of differential fi — b cross sections, where the muon is from a 
semi-leptonic b decay and the b is identified using precision track reconstruction in jets. 
The semi-differential correlated cross sections, dcr/dE T , dcx/dp T , and d<r/d<50(/i — b) for 
> 9 GeV/c, |?7 M | < 0.6, E T > 10 GeV, \rf \ < 1.5, are presented and compared to 
next-to-leading order QCD calculations. 

PACS Numbers: 13.85. Qk, 13.87.-a, 14.65.Fy 

1 Introduction 

Measurements of b production in pp collisions provide quantitative test of perturbative QCD. 
Single integral b cross section measurements at y/s = 1.8 TeV have been systematically 
higher than predictions from next-to-leading order (NLO) QCD calculations ]I|, 0, [J. These 
cross section measurements, from inclusive b — > lepton decays and exclusive B meson decays 
(B + ^J/ipK + ), use the kinematical relationship between the decay product (e.g, the lepton) 
and the b quark spectra to obtain the production cross section integrated over a rapidity 
range | y | < 1 and a range from a threshold p™ n to infinity. Single differential B meson 
cross section measurements Q are also systematically higher than the NLO prediction. 

Semi-differential b — b cross sections give further information on the underlying QCD 
production mechanisms by exploring the kinematical correlations between the two b quarks. 
Comparison of NLO predictions with experimental measurements can give information on 
whether higher order corrections serve as a scale factor to the NLO prediction or change the 
production distributions. As future high precision B decay measurements at hadron colliders 
{e.g., CP violation studies in B° ^J/ipK® 0) may depend upon efficient identification of the 
decay products of both b quarks, understanding of the correlated cross sections is necessary. 

This paper describes measurements of /i — b correlated cross sections as a function of the 
jet transverse energy (da/dEx, where Ex = E x sinO) and transverse momentum (da/dpx) 
of the b and as a function of the azimuthal separation (dcr/d<50) between the muon and b 
jet, for p^ > 9 GeV/c, |^| < 0.6, E T > 10 GeV/c, \rf \ < 1.5. The data are 15.08 ± 0.54 
pb _1 of pp collisions at yfs = 1.8 TeV collected with the CDF detector between August, 
1992 and May, 1993. We make use of two features of B hadrons to separate them from the 
large jet backgrounds at 1.8 TeV: the high branching fraction into muons (~ 10% |8[]) and 
the relatively long lifetime (~ 1.5 picoseconds ||). The advent of precision silicon microstrip 
detectors, with hit resolutions approaching 15 /zm, provides the ability to efficiently identify 
the hadronic decays of B hadrons as well as the semi-leptonic decays. 

We use the identification of a high transverse momentum muon as the initial signature 
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of the presence of b quarks. In pp collisions, high transverse momentum muons come from 
the production and decay of heavy quarks (c, b, t), vector bosons (W,Z°), and light mesons 
(ti,K). Additional identification techniques are necessary to convert a //— jet cross section 
into a fi — b cross section. 

For these measurements, the first b is identified from a semi-leptonic decay muon and 
the other b (referred to for simplicity as the 6, though we do not perform explicit flavor 
identification for either b) is identified by using precision track reconstruction in jets to 
measure the displaced particles from b decay. Jets are identified as clusters of energy in the 
calorimeter ||. In this paper, a jet energy (or jet transverse energy) refers to the measured 
energy in the cluster. A procedure to simultaneously unfold the effects of detector response 
and resolution is used to translate the results from b jets to b quarks. 

It should be noted that we have chosen to report the measurements as differential fi — b 
cross sections rather than b — b cross sections in order to facilitate comparison to calculations 
of the production cross sections. The process of converting a muon cross section to a quark 
cross section includes systematic uncertainties [0] with strong dependence on both produc- 
tion, fragmentation, and decay models. By presenting fi — b cross sections, we facilitate the 
future comparison of the experimental results to different models, since the data results and 
uncertainties are not tied to specific models. 

Section ^| describes the detector systems used for muon and b jet identification. Section |3] 
contains descriptions of the muon and jet identification requirements. The b jet counting 
is discussed in section f|. In section [5], the muon and b jet identification efficiencies and 
acceptances are described. The cross section results, the calculation of additional physics 
backgrounds, and jet to quark unfolding are discussed in section ||. Section [7| closes with a 
discussion of the experimental results. 

2 Detector Description 

The CDF has been described in detail elsewhere [|Kj. The analysis presented in this paper 
depends on the tracking and muon systems for triggering and selection, while identification 
of hadronic jets uses the information from the calorimeter elements. 

2.1 Tracking and Muon Systems 

This analysis uses the silicon vertex detector (SVX) ||11|| , the vertex drift chamber (VTX) 
and the central tracking chamber (CTC) [|12| for charged particle tracking. These are all 
located in a 1.4 T solenoidal magnetic field. The SVX consists of 4 layers of silicon-strip 
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detectors with r — readout, including pulse height information, with a total active length 
of 51 cm in the range —27.3 < z < 27.3 cm 0. The pitch between readout strips is 
60 /im on the inner 3 layers and 55 /im on the outermost layer. A single point spatial 
resolution of 13 /im has been obtained. The first measurement plane is located 2.9 cm from 
the interaction point, leading to an impact parameter resolution of ~ 15 /im for tracks with 
transverse momentum, pp, greater than 5 GeV/c. The VTX is a time projection chamber 
providing information out to a radius of 22 cm and | 77 |< 3.5. The VTX is used to measure 
the pp interaction vertex (zq) along the z axis with a resolution of 1 mm. The CTC is a 
cylindrical drift chamber containing 84 layers, which are grouped into alternating axial and 
stereo superlayers containing 12 and 6 wires respectively, covering the radial range from 28 
cm to 132 cm. The momentum resolution of the CTC is Spt/pt = 0.002 x for isolated 
tracks (where is in GeV/c). For tracks found in both the SVX and CTC, the momentum 
resolution improves to 0.0009 x p T ® 0.0066 (where is in GeV/c). 

The muon system consists of two detector elements. The Central Muon system (CMU) [|T3]], 
which consists of four layers of limited streamer chambers located at a radius of 384 cm, be- 
hind ~ 5 absorption lengths of material, provides muon identification for the pseudorapidity 
range \r)\ <0.6. This 77 region is further instrumented by the Central Muon upgrade system 
(CMP) ||14|| , which is a set of four chambers located after « 8 absorption lengths of mate- 
rial. Approximately 84% of the solid angle of \rj\ <0.6 is covered by the CMU, 63% by the 
CMP, and 53% by both. Muon transverse momentum is measured with the charged tracking 
systems and has the tracking resolutions described above. CMU (and CMP) segments are 
defined as a set of 2 or more hits along radially aligned wires. 

2.2 Calorimeter Systems 

This analysis uses the CDF central and plug calorimeters, which are segmented into separate 
electromagnetic and hadronic compartments. In all cases, the absorber in the electromagnetic 
compartment is lead, and in the hadronic compartment, iron. The central region subtends the 
range \r]\ < 1.1 and spans 2ir in azimuthal coverage, with scintillator as the active medium. 
The plug region subtends the range 1.1 < |r?| < 2.4 with gas proportional chambers as the 
active media, again with 2ir azimuthal coverage. The calorimeters have resolutions that 
range from 13.7%/v / E T " © 2% for the central electromagnetic to 106%/v / Er © 6% for the 
plug hadronic fTSfl . 
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2.3 Trigger System 



CDF uses a three-level trigger system [ 16| . Each level is a logical OR of a number of triggers 
designed to select events with electrons, muons or jets. The analysis presented in this paper 
uses only the muon trigger path. Section ^ includes a description of the trigger efficiencies 
for muons. 

The Level 1 central muon trigger requires a pair of hits on radially aligned wires in the 
CMU system. The p^ of the track segment is measured using the arrival times of the drift 
electrons at the wires to determine the deflection angle due to the magnetic field. The 
trigger requires that the segment have pt > 6 GeV/c, with at least two confirming hits in 
the projecting CMP chambers. 

The Level 2 trigger includes information from a list of r — tracks found by the central 
fast tracker (CFT) []17 |, a hardware track processor which uses fast timing information from 
the CTC as input. The CFT momentum resolution is Sp^/pT ~ 0.035 x p^, with a plateau 
efficiency of 91.3±0.3% for tracks with above 12 GeV/c. The CMU chamber segment is 
required to match a CFT track with > 9.2 GeV/c within 5° in the </> coordinate. 

The Level 3 trigger makes use of a slightly modified version of the offline software recon- 
struction algorithms, including full 3 dimensional track reconstruction. The CMU segment 
is required to match a CTC track with p^ > 7.5 GeV/c, extrapolated to the chamber radius, 
within 10 cm in r — 0. Confirming CMP hits are required. 



3 Dataset Selection 

Beginning with the sample of muon triggered events, we select events with both a well 
identified muon candidate and a minimum transverse energy jet. A primary vertex is found 
by a weighted fit of the VTX zq vertex position and SVX tracks. An iterative search removes 
tracks with large impact parameters (the distance of closest approach in the r — <fi plane) from 
the fit. Since the b jet identification technique (described in section f|) depends upon the 
precision track reconstruction in the SVX, we require the event primary vertex | zq \ < 30 cm. 
In this section, we discuss the identification variables, efficiency, and geometric acceptance 
for muon and b jet candidates. Table ^| contains a summary of the muon efficiency and 
acceptance results and table f| contains a summary of the b jet identification and acceptance 
results. 



7 



3.1 Muon Identification 



Muons are identified as a well matched coincidence between a track in the CTC and segments 
in both the CMU and CMP muon systems. The CTC track is required to have > 9 
GeV/c and point back to within 5 cm in z of the found primary vertex. The measured 
track is extrapolated to the muon chambers and is required to match the muon chamber 
track segment position to < 3a in the transverse direction (for both CMU and CMP) and 
< \J~Y2o in the longitudinal direction (for CMU). In all cases, a includes the contributions 
from smearing due to multiple scattering in the absorber and the muon chamber resolution. 
We require that the track be found in the SVX. 

There are 144097 events passing all the muon requirements in this data sample. In the 
case where there is more than 1 identified muon in an event, we take the highest pt muon as 
the b candidate muon. The fraction of muons from b decay is measured to be approximately 
40% [15], with a fraction from charm decays of approximately 20%. Figure [I] shows the 



transverse momentum spectrum for the muons in this dataset. The flattening of the slope 
at high is due to muons from electroweak boson decay. 

3.2 Jet Identification 

Jets are identified in the CDF calorimeter systems using a fixed cone (in rj — <fi space) 
algorithm. A detailed description of the algorithm can be found in reference @. For this 
analysis, we use a cone radius of 0.4. We require that jets have transverse energy, Et = E x 
sin 9 (where E is the total energy in the cone), greater than 10 GeV, and | rj |< 1.5. There 
are 50154 events passing the muon and jet Et requirements. We use tracking techniques to 
identify b jets, so the pseudorapidity range is restricted to the region with tracking coverage. 
All jet energies in this paper are measured energies, not including corrections for known 
detector effects(e.#., calorimeter non-linearities). An unsmearing procedure, described in 
section is used to convert measured jet Et distributions to parton momentum distributions. 

We associate SVX tracks to a jet by requiring that the track be within the cone of 0.4 
around the jet axis. To remove tracks consistent with photon conversions and K$ or A 
decays originating from the primary vertex, we require that the impact parameter, d, be less 
than 0.15 cm. In addition, track pairs consistent with Kg — > tt + tt~ or A — ► pit decays are 



removed. We select jets with two or more well measured tracks pr > 1 GeV/c, with 
positive impact parameters. The impact parameter sign is defined to be +1 for tracks where 
the point of closest approach to the primary vertex lies in the same hemisphere as the jet 
direction, and —1 otherwise, for 



S 




Figure 1: The /i p T spectrum for the 9 GeV/c sample. There are 144097 events, with 80 
having p T > 60 GeV/c. The enhancement above 25 GeV/c is due to the presence of muons 
from W and Z boson decays. 
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We require that the distance, AR, in rj — <p space between the muon and the jet axis be 
greater than 1.0. There are 16842 events passing all the muon and jet requirements. The 
AR separation is chosen so that tracks clustered around the jet axis are separated from the 
muon direction, in order to have physical separation of the b and b decay products. As there 
may be more than one jet in an event passing these requirements, we select the jet with the 
lowest jet probability (defined in section so as to have a unique combination of /i — jet 
in each event. 



4 b Jet Counting 

The b jet is not identified on an event-by-event basis, but instead by fitting for the number 
of b jets present in the sample. For each jet, we combine the impact parameter information 
for tracks in the jet cone into one number which describes the probability that the given 
collection of tracks has no decay products from long lived particles. In a b jet, there will be 
a significant number of tracks from the B hadron decay, and hence the probability for a b 
jet will be much less than 1. 



4.1 The Jet Probability Algorithm 



The b jet identification makes use of a probability algorithm [|18j which compares track impact 
parameters to measured resolution functions in order to calculate for each jet a probability 
that there are no long lived particles in the jet cone. This probability is uniformly distributed 
for light quark or gluon jets (we refer to these jets as prompt jets), but is very low for jets 
with displaced vertices from heavy flavor decay. We now briefly describe the transformation 
from the track impact parameters to the jet probability measure. 

The track impact parameter significance is defined as the value of the impact parameter 
divided by the uncertainty in that quantity, which includes both the measured uncertainties 
from the track and primary vertex reconstruction. Figure ^| shows the distribution of impact 
parameter significance (sq = d/cr) from a sample of jets taken with a 50 GeV jet trigger [[TBI , 



overlayed with a fitted function. The tails of the distribution come from a combination of 
non-Gaussian effects and true long lived particles. Using a combination of data and Monte 
Carlo simulation of heavy flavor decays, we estimate approximately 30% of the tracks with 
| so |> 3.0 are from the decay products of long lived particles, which is consistent with the 
excess in the positive s side of the distribution. The negative side of the fitted function, 
R(s), is used to map the impact parameter significance Sq to a track probability measure: 
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I-oo R ( S ) ds 

The track probability is a measure of the probability of getting a track with impact parameter 
significance greater than s. The function R(s) can be defined for both Monte Carlo simulated 
datasets and the jet dataset. The mapping of the resolution function to the track probability 
distribution removes differences in the resolution between the simulated detector performance 
and the true detector performance and creates a variable which is consistent between the 
two datasets. 

The jet probability is then calculated from the independent track probabilities as: 

= n E LJ ^ HL > (2) 

where 

n = p 1 p 2 ■ ■ ■ p N (3) 

is the product of the individual probabilities of the selected tracks. For the rest of this paper, 
when the track selection requirements pick tracks with negative signed impact parameters, 
we will refer to the measure as the "negative jet probability". When the track selection 
requirements pick tracks with positive signed impact parameters, we will refer to the measure 

aS Pjet- 

Figure |3| shows the distribution of the negative jet probability in the 50 GeV sample. Since 
this distribution reflects the smearing of the impact parameter significance distribution due 
to resolution effects, we expect that the Pj et distribution for prompt (light quark and gluon) 
jets to be similar. Simulated jets containing heavy flavor decays show distinct differences 
from this distribution, peaking at low values of Pjet- in figure [|, we show the distributions 
of logio {Pjet) for b, charm, and prompt jets. 

We have found that the Pj et shape for heavy flavor jets is affected by the number of 
tracks used in the calculation of Pj et which are also used in the primary vertex fit. The 
turnover visible in the b and charm distributions around -3 in logio (Pjet) is a combination of 
the vertex requirements [dja < 3 for tracks in the fit) and the b and charm lifetimes, b and 
charm jets are affected differently, due to differences in lifetime and decay multiplicities. 

4.2 b Jet Fit Technique 

We use a binned maximum likelihood fit to distinguish the b, c, and prompt jet contributions 
in the sample. For a binned likelihood fit, we find that logio (-Pjet) shows stronger differen- 
tiation between b, c, and prompt jets (see figure £|) than Pj et and use this variable in the 
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Figure 2: A sample track resolution function, including fits to both positive and negative 
signed impact parameters. The function is fit to 2 gaussians plus two exponentials, one for 
the positive side and one for the negative side. The excess on the positive side is attributable 
to long lived particles in the sample. 
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Figure 3: The negative jet probability spectrum, calculated using tracks with negative signed 
impact parameters, in a sample of 50 GeV jets. 



13 



fitting algorithm. We fit over the range -10 — in logio(-Pjet), where the b, c, and prompt 
contributions are constrained to be positive. No other constraints are included in the fit. 
We model the prompt jets with an exponential distribution, since a logarithm transforms a 
uniform distribution to an exponential distribution. 

We have explored the effect of different Monte Carlo samples to construct the input shape 
used in the fit. Using different input b jet Monte Carlo samples (see section |5|) compared to 
the test distribution shows a 5% change in the fit fractions. Changing the average b lifetime 
by 6% ]19j changed the fit fraction by 3%. We include a 5.8% systematic uncertainty to 



our fit results to account for systematic uncertainties in the fitting procedure and uncertainty 
on the b lifetime. 

In figure || we show the distribution of logio(P Je t) for all jets, E T > 10 GeV, in the muon 
sample, overlayed with the fit results. In this sample, the fit finds 2484 ± 94 b jets, 1988 
± 175 c jets, and 12368 ± 157 prompt jets for a total of 16840. There are 16842 events in 
the data sample. Figure |6] shows three comparisons of the data and fit results, showing the 
bin-by-bin difference in the results, the bin-by-bin difference divided by the errors, and the 
distribution of the difference divided by the errors. In these distributions, the errors are the 
statistical errors in the data points. We do not include any error on the Monte Carlo shapes. 
From these distributions, we can see that the inputs model the data well. The difference 
divided by the errors has a mean of 0.04 and RMS of 0.95. 

For the semi-differential measurements, we do an independent fit of the logio(-Pjet) distri- 
bution and then correct for the acceptance in each Et or 5<j) bin. Table [3] contains a summary 
of the number of total jets and the number of b jets in each Et and 5(f) bin considered. 

5 Acceptance and Efficiency 
5.1 Muon Requirements 

The muon geometric acceptance is the fraction of events with a muon in the good fiducial 
region of the CMU and CMP chambers, starting from a sample where the muon has > 
9 GeV/c and | rj |< 0.6. Note that this term is only a geometric acceptance and does not 
include kinematical cuts on the muon. 

The geometric acceptance is studied with a b — > /i Monte Carlo generator (which in- 
cludes the sequential decays b — > cfi), with the input spectra coming from the next to 
leading order calculation of b — b production by Mangano, Nason, and Ridolfi (MNR) [^2| . 



The input spectra use the MRSD0 structure functions p3f and renormalization scale /zq 



'ml + (p^ +p T 2 2)/2, with = 4.75 GeV/c 2 . This generator produces b quarks and B 
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Figure 4: The logio(-Fjet) distributions used as inputs to the fitting program. The b and c 
shapes are smoothed versions of Monte Carlo distributions, while the primary shape is an 
exponential function. The three distributions are normalized to equal area and shown on 
the same vertical scale. 
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Figure 5: For all jets (Ex > 10 GeV) in the fi sample, we show the data distribution overlayed 
with the fit results. There are two events in the data with logio(Pjet) < -10. Statistical errors 
on the data and the fit results are included. The fit results model the data well over the 
entire range of the fit. 
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Figure 6: Various comparisons of the data distribution and the fit results. We show (a) the 
bin by bin difference between the data and the fit results, (b) the bin-by-bin difference scaled 
to the errors, and (c) the distribution of the difference scaled to the errors, with mean 0.04 
and RMS = 0.95. In all cases, the errors are the statistical error in the data points and the 
fitted results. 
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Et Range 


Number of Jets 


Estimated Number 






of b Jets 


10 — 15 


5174 


547 ± 49 


15 — 20 


3818 


618 ± 47 


20 — 25 


2553 


a r o i on 

453 ± 39 


o cr on 

25 — 30 


1698 


oto i on 

278 ± 30 


oil — 4U 




Q07 -U QQ 
OZ( ± oo 


40 — 50 


819 


1 y|n I on 

140 ± 20 


rn 1 n r\ 

50 — 100 


o A n 

849 


1 nT i 1 n 

107 ± 19 


50 Range 






— - 

w 8 


43 


4.8 +5.5-4.8 


7T 7T 

8 4 


83 


25.0 ± 8.6 


7T 07T 
4 8 


230 


54.7 ± 13.3 


37T 7T 


336 


78.2 ± 15.9 


8 . 2 

7T 07T 

2 8 


519 


105. ± 18.5 


57T 37T 

8 4 


1008 


160. ± 25. 


37T 77T 

1 8 


3229 


461. ± 42. 


^ — 7T 
8 


11394 


1593. ± 75. 



Table 1: b fit results as a function of jet Et and S(p between the muon and b jet. We have 
not included a common systematic uncertainty of 5.8%. 



hadrons, using the Peterson fragmentation form with e = 0.006 ± 0.002 ||25| . B 
hadrons are decayed according to the CLEO Monte Carlo program, QQ pB| . We select 
events with a6->/j decay, with muon > 9 GeV/c and | 77 |< 0.6. 

For these studies, event vertices are distributed along the Gaussian with mean 



= -1.4 cm and a = 26.65 cm ||21|| , which is a good approximation to the average conditions 
seen in the data. The muons are propagated to the CMU and CMP chamber radii, including 
the effects of the central magnetic field and multiple scattering. The acceptance is then 
defined as the fraction of muons which are in the good fiducial area of both the CMU and 
CMP chambers and is found to be 53.0 ± 0.3% (statistical), independent of variations of the 
e parameter from 0.004 to 0.008. 

The muon trigger and selection depends significantly upon the track reconstruction ef- 
ficiency in the CTC. We have defined our efficiencies to be multiplicative, so that we can 
measure them independently. In this section, the efficiencies of the individual selection re- 
quirements, and methods of measuring them, are described. 

The trigger efficiency is measured using independently triggered samples for each level of 
the system, where the efficiency is expressed as a function of the muon p^- Figure [7] shows 
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the efficiency curves for the 3 levels of the trigger system. The efficiency curves are then 
convoluted with the px spectrum of the muons, to extract the efficiency for a muon with 
Pt > 9 GeV/c. This convolution is done independently for the differential Et cross section 
bins (see table 0), since the muon p^ spectrum may depend upon the transverse momentum 
distribution of the 6 jet recoiling against the b — > \i decay. For b jets with Et > 10 GeV, the 
combined LI, L2, and L3 trigger efficiency is measured to be 83.0 ± 1.7 %. 

The vertex requirement, | z |< 30 cm, is studied in a minimum bias trigger dataset, 
comparing the vertex distribution to the predicted shape, including the measured longitudi- 
nal distribution of the proton and anti-proton bunches and the effects of the accelerator (3 



function ||21|| . The efficiency is found to be 74.2 ± 2.1 %, where the uncertainty comes from 
uncertainty in the measured beam longitudinal distributions and f3 function. 

The track finding efficiency in the CTC is a function of the density of charged particles. By 
embedding Monte Carlo simulated track hits into data samples, we quantify the probability 
of finding the Monte Carlo simulated track as a function of the relative density of CTC 
hits. The quantified probability is convoluted with the hit density distribution for the muon 
sample. The track finding efficiency is measured to be 96 ± 1.7 %, where the uncertainty 
represents the change in the result using different parametrizations of the probability curve 
vs hit density. 

The combined x 2 matching efficiency is measured in a J/ip —>■ /i + /U~ sample identified 
by tracking and mass requirements and is found to be 98.7 ± 0.2 %. The muon segment 
reconstruction efficiency is found to 98.1 ±0.3 %, resulting in a combined efficiency of 96.8 ± 
0.4 %. 

The track finding efficiency in the SVX is studied in the 9 GeV/c muon sample, requiring 
the CTC track to extrapolate to a good SVX fiducial region. The efficiency is found to be 
90 ± 1%, where the uncertainty is the statistical error only. 

5.2 b Jet Requirements 

The b jet acceptance combines the fiducial acceptance of the SVX and the CTC, the track re- 
construction efficiency, and fragmentation effects and the AR separation requirement. These 
tracking and AR effects are studied separately, with a full simulation used for the combina- 
tion of the track requirements and fiducial acceptance, while a MNR based \l — b model is 
used for the AR acceptance. The b jet acceptance is calculated separately as a function of 
the jet Et and azimuthal opening angle between the muon and the jet. 

Monte Carlo samples for b and c quarks are produced using ISA JET version 6.43 |[27|| . 
The CLEO Monte Carlo program |26| is used to model the decay of B hadrons. b quarks 
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Figure 7: The trigger efficiency curves for the 3 levels of the trigger system. The trig 
efficiency is the product of the three curves, convoluted with the \x pp spectrum. 
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Et bin 


Trigger Efficiency 


10 — 15 GeV 
15 — 20 GeV 
20 — 25 GeV 
25 — 30 GeV 
30 — 40 GeV 
40 — 50 GeV 
50 — 100 GeV 


82.6 ± 1.7 % 
83.0 ± 1.7 % 
83.4 ± 1.7 % 

83.6 ± 1.7 % 

83.8 ± 1.7 % 

83.9 ± 1.7 % 

83.7 ± 1.7 % 


All E T 


83.0 ± 1.7 % 



Table 2: The muon trigger efficiency for each jet Et bin. A common 2% uncertainty is 
assigned to each bin. 



Geometric Acceptance 


53.0 ± 0.3 % 


CTC Track Finding 
Matching Efficiency 
Z Vertex Requirements 
SVX Track Finding 


96.0 ± 1.7 % 
96.8 ± 0.4 % 
74.2 ± 2.1 % 
90 ± 1 % 


Combined Acceptance 
and Efficiency 


32.9 ± 1.1 % 



Table 3: Summary of muon acceptance and efficiency numbers. The trigger efficiency is 
applied on a bin by bin basis for the jet Et measurement. 



produced using the HERWIG Monte Carlo |23j and PYTHIA Monte Carlo J^J programs 
are also used for systematic studies. The ISAJET and PYTHIA samples used the Peterson 
form as the fragmentation model, with e = 0.006 ± 0.002. While none of these generators 
use a NLO calculation of b production, the r\ distribution of the quarks agrees well with the 
NLO calculation. For tracking efficiency studies, events with a muon with > 8 GeV are 
passed through the full CDF simulation and reconstruction package. The simulation used 



an average b lifetime of cr = 420 /im [20 



The track acceptance represents the fraction of b quarks, Et > 10 GeV, \q\ < 1.5 which 
produce jets with at least 2 good tracks inside a cone of 0.4 around the jet axis, where 
there is also a b quark which decays to a muon with pt > 9 GeV within the CMU-CMP 
acceptance. The average track acceptance for the b is 51.4 ± 0.8%. It ranges from 45.7 ± 
1.1% (statistical error only) for 10 < E T < 15 GeV to 64.8 ± 2.6% for 50 < E T < 100 GeV. 

We have compared the values for the b track acceptance from ISAJET samples to the 
acceptance from HERWIG samples. The acceptance agrees within the statistical error in 
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Et Range 


Track Acceptance 


AR Acceptance 


10 — 15 
15 — 20 
20 — 25 
25 — 30 
30 — 40 
40 — 50 

ou — 1UU 


45.7 ± 1.1 ± 2.3 % 
55.9 ± 1.7 ± 2.8 % 
58.1 ± 2.5 ± 2.9 % 
61.3 ± 3.5 ± 3.1 % 

61.7 ± 3.8 ± 3.1 % 

64.8 ± 2.6 ± 3.2 % 

00. U IE Z.u IE O.O /o 


86.9 ± 1.0 m % 

88.2 ± 1.5 til % 

oo o i o r\ -f-2.2 ny 

88.3 ± 2.0 I 2 .o % 
88.3 ± 2.3 tti % 
87.9 ± 3.4 tH % 
87.1 ± 3.5 tti % 

ok e _i_ q 7 +5.2 (w 
OO.O IE o. ( _i g /o 


5(f) Range (radians) 






— - 

7T 7T 

8 4 

7T on 
4 8 

8 " 


46.3 ± 1.4 ± 2.6 % 

47.3 ± 1.4 ± 2.6 % 

51.4 ± 0.8 ± 2.6 % 
51.4 ± 0.8 ± 2.6 % 


6.9 ± 0.03 % 
20.8 ± 0.2 til % 
74.7 ± 0.9 til % 
100 % 



Table 4: fi — b track and AR acceptance as a function of jet Ex and 5(f) (statistical and 
systematic uncertainties). There is a common (relative) systematic uncertainty of 5% in the 
tracking efficiency For 5(f) > 1 radian, the AR acceptance is 100% by definition. 

the samples as a function of Ex, differing at the 5% level. We include this variation as 
an additional systematic uncertainty on the track acceptance. Comparisons of inclusive jet 
track acceptances from an IS A JET sample and from data show reasonable agreement. 

For the calculation of the AR acceptance, we have used a model based on the MNR 
calculation fl22" |. This calculation can be used to give exact 0(a^) results in situations where 
kinematical cuts have been applied at the parton level. We have made additions to the 
calculation to model the /x — b differential cross sections. 



The MNR calculation produces the vectors p b , p b , and p9 luon with appropriate weights. 
We include additional weighting for the following: 

• Probability of > 9 GeV/c for given p T , "P(pt)Pt) 

• Probability of E^ jet in a given Et bin for given p T , P(E T ,p T ) 

^(PtjPt) i s defined as the fraction of b quarks, with given p T , which decay into muons 
with > 9 GeV/c. We use the b — > /i Monte Carlo generator described above to derive this 
function, using B(b — > fi) = 0.103 ± 0.005 |§. Since the probability is defined as a function of 
p T , the exact shape of the p\ distribution does not enter into the result. Figure || shows the 
value P(pt)Pt) as a function of p\. The three curves are for different values of the Peterson 
e parameter used in the fragmentation model. In addition to this probability weighting, we 
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also smear the b quark direction in pseudo-rapidity and azimuth. The smearing is based on 
the results from the b — *■ \i Monte Carlo generator. 

V(E^,p\-) is defined as the probability that a b quark, with given p T , would produce a 
jet with given E T . Using the methods outlined in section |6], we have a binned probabil- 
ity distribution in Et for each p\. Since the measured jet Et integrates over a range in 
pseudo-rapidity and azimuth (a cone of radius 0.4), we approximate this clustering effect by 
clustering partons (adding the b and the gluon momenta vectorially) within the same cone 
size. For the rest of this paper, when we discuss the E T or p T theory distributions, it means 
the clustered partons E T or p T . 



We use a renormalization and factorization scale /i = ymf + (p^ + p^)/2, MRSA struc- 
ture functions ||30|| , and mb = 4.75 GeV/c 2 . Applying the additional weights and the appro- 
priate kinematical cuts (| rf |< 0.6 and | rf |< 1.5), we obtain the calculated d<7/dE T and 
da/d5(p(fi — b) distributions. We create the same distributions with the requirement that 
the muon and b be separated by AR > 1 and do a bin by bin comparison of the calculated 
cross sections to define the AR acceptance. We have varied the renormalization scale, b 
quark mass, and parton distribution functions used in the MNR calculation to estimate the 
systematic uncertainties in the AR acceptance. Table £| shows the bin by bin values used in 
the differential cross section measurements. 



6 Cross Section Results 

The cross section results are presented as fi — b cross sections. Since we have not specifically 
done flavor identification, there is an additional factor of 1/2 in the calculation of the cross 
sections. For the semi-differential measurements, we do an independent fit of the logio(-Pjei) 
distribution and then correct for the acceptance in each Et or 5<f) bin. With the number 
of b jets from table [l], the bin by bin trigger efficiencies from table ^ the combined muon 
acceptance and efficiency from table |3|, and the b track and AR acceptances from table [|, 
we calculate the the cross section in each Et and 8cj) bin considered. The sum of the 7 Et 
bins is 614.4 ± 63.0 pb and the sum of the 8 6(f) bins is 633.0 ± 70.6 pb. The results are 
summarized in table |5|. 

6.1 Physics Backgrounds 

There are backgrounds which need to be included before comparing to theoretical predictions 
on b — b production, since there are additional sources of pi — b production. Specifically, the 
decay products of light mesons (n, K) produced in association with b — b pairs or heavy 
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Figure 8: The probability of a b — > /j, decay, with pj, > 9 GeV/c, as a function of p^- We 
have included the branching fraction B (b —>//) = 0.103. The curves represent three choices 
of the Peterson e parameter used in the fragmentation process. 
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Et 


Range 


Cross Section (pb) 


10 


- 15 


168.1™ 


15 


- 20 


152.2_ 12 ; 8 


20 


- 25 


106.6I 10 . 5 


25 


- 30 


61.93+^ 


30 
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79 cq+8.96 


40 


- 50 


29.81+1™ 


50 


- 100 


r,o -i q+4.38 
ZO.lO_4.23 


5(f) Range 






7T 

8 

7T 
4 

3?r 
8 


7T 

8 

7T 

4 

3tt 
8 

7T 

2 


30.86±Jg:S 
17.30+|;|? 

18.48 ± 3.76 


TV 

2 


5tt 


24.81 ± 4.37 


5ir 
8 


3^7T 

4 


37.81 ± 5.91 


3n 
4 


7tt 
8 


108.9 ± 9.93 


7tt 

8 


7T 


376.4 ± 17.7 



Table 5: /i — b cross sections as a function of jet Ex and S(j) between the muon and b jet. 
We have not included a common systematic uncertainty of 9.3% in the results. Physics 
backgrounds have not been subtracted at this stage. 



particles (e.g, the Z° boson, top quark production) can give a similar signature. 

A contribution to the sample occurs when the identified muon is not coming from a b 
quark decay but instead from the decay of a light meson (ir or K) or charm quark. In 
the inclusive muon sample, the b fraction is measured to be approximately 40% [|15|], with a 
charm fraction of approximately 20% and the remaining 40% from the decay of light mesons. 
Since jets from gluons are the dominant production process in this jet Et range, we assume 
that the light mesons come predominantly from gluon jets. With the further assumption 
that the gluon splitting to bb probability is approximately 1.5% p|, we estimate that in 
0.6% (0.015 x 0.4) of the muon events we correctly identify the b but the muon is from a 
light meson decay. The case where the identified muon comes from the decay of a charm 
particle can be estimated in a similar manner. With the same assumptions about the gluon 
splitting to heavy quark probability (1.5%), a measured charm fraction of 20%, and that 
approximately 75% of charm quarks are produced via gluon splitting, we estimate that in 
0.2% (0.015 x 0.75 x 0.2) of the muon events we correctly identify the b but the muon is from 
a charm particle decay. 

With an identified fraction of 40% b muons and 50% of the produced 6's from gluon 
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Et Range 


Cross Section (pb) 




Statistical Uncertainty only 


10 — 15 


0.43 ± 0.06 


15 — 20 


0.75 ± 0.08 


20 — 25 


r\ o c\ i r\ r\r\ 
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— ou 


n fin -i- n n7 

U.OU it U.U ( 
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ou — 4tU 


n S7 4- n m 
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40 — 50 


0.12 ± 0.02 


50 — 100 


0.015 ± 0.008 


5(f) Range 




— - 

w 8 





7T 7T 

8 4 





7T 37T 

4 8 


0.015 ± 0.014 


3-7T 7T 


0.031 ± 0.012 


8 S 2 
7T 07T 

2 8 


0.036 ± 0.013 


8 4 


0.11 ± 0.024 


i 8 


0.53 ± 0.05 


8 " 


2.88 ± 0.12 



Table 6: Contributions from Z° — > fib to the cross section as a function jet Et and 8(f) 
between the ji and b jet. There is an addition 8.0% uncertainty in the overall normalization. 



splitting [3"3~|, in 20% of the muon events we correctly identify the b and the muon from the b 
decay. Combining these calculations yields a fractional background in the \x — b cross section 
of 0.04 (= (0.006 + 0.002)/0.20). We assume that this background has the same shape as 
the signal and reduce the cross sections by a constant 4.0 ± 2.0% (the uncertainty is taken 
as half the change) . 

We have used the PYTHIA Monte Carlo program to generate Z° — > bb events, and the 
CLEO Monte Carlo program for the decay of the resulting B hadrons. We normalize the 
production cross section to measured CDF cross section of Z° — > e + e~ M, RTJ , and apply 
the same fi and b jet requirements as presented in section [| The predicted cross section 
remaining after these requirements is 3.6 ± 0.28 pb, where the uncertainty includes the 
relative normalization to the dielectron decay mode, the b — > fi branching fraction, and 
acceptance uncertainties. Table [| shows the contributions from this process in the same Et 
and S(j) bins as in table ||. 

Top quark production and decay can also contribute to the \i — b cross sections. The 
CDF measurement of the total top cross section is 6-8-2! pb [P^] . However, once we account 
for branching fractions and acceptance criteria, the total cross section from this process is 
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less than 1 pb and will not be considered further. 
6.2 Jet Unsmearing Procedure 

The cross sections measured above depend upon the selection of jets with E T > 10 GeV, and 
in the case of the dcx/dEx distribution, depend upon the binning of the distribution. Jets 
coming from b quarks with transverse momentum p\ will contribute to more than one bin 
in the measured distribution, due to the combined effects of calorimeter energy response, 
calorimeter energy resolution, and quark fragmentation. An unsmearing procedure has been 
developed at CDF to account for these effects. 

We use Monte Carlo produced samples to define the expected jet Et response distri- 
bution for a given quark pp. An iterative procedure is used to correct the measured cross 
sections. The quark distribution is described by a smooth function and smeared with 
the simulation derived E T response functions. The input distribution is adjusted until the 
smeared distribution matches the measured distribution. We then perform a simultaneous 
unfolding of the measured jet Et spectrum to the parton pr spectrum to account for energy 
loss and resolution. This unfolding corrects both the cross section and Et (pr) axes. 

6.2.1 Response Functions 

The calorimeter single particle response in the range 0.5 to 227 GeV has been determined 
from both test beam data and isolated tracks from collider data. A Monte Carlo simulation 
incorporating the calorimeter response and the ISAJET, HERWIG, and PYTHIA samples 
is used to determine a response function for b jets in the Et range 5 to 150 GeV, including 
energy loss, resolution, and jet finding efficiency effects. For each pr, the response function 
represents the probability distribution for measuring a particular value of Et- These response 
functions are convoluted with the expected b pp distributions, creating an expected Et 
distribution. 

6.2.2 Unsmearing 

The input b distribution comes from the // — b model described in section [|, where we 
have required a muon with pr > 9 GeV/c. We have parametrized the distribution with a 
multi-quadric function and varied a scale parameter until the smeared distribution matches 
the measured distribution. Figure ^| shows the best match b distribution, overlayed 
with the smeared distribution. Table [7] shows the unfolding effects on the cross section and 
transverse momentum. Note that the unsmearing procedure introduces correlated systematic 
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Jet E T Bin 


Mean Jet Ex 


a (pb/GeV) 


Mean b px 


a (pb/GeV/c) 


10 — 15 


12.38 


32.20 


25.28 


27.66 


15 — 20 


17.35 


29.07 


30.67 


24.62 


20 — 25 


22.30 


20.30 


35.99 


18.78 


25 — 30 


27.34 


11.77 


41.20 


12.52 


30 — 40 


34.31 


6.88 


48.38 


7.13 


40 — 50 


44.36 


2.85 


59.00 


3.05 


50 — 100 


63.19 


0.44 


79.18 


0.57 



Table 7: Smeared and unsmeared means and cross sections for the 7 bins in the differential 
Px measurement. The cross sections are after background subraction and are presented here 
without uncertainties. Note that the unsmearing procedure introduces correlated uncertain- 
ties in the bins. 

uncertainties in the bins. 

6.2.3 Systematic Uncertainties 

Systematic uncertainties in the smearing procedure arise from uncertainties in the knowledge 
of the calorimeter energy scale, the calorimeter resolution, the jet finding efficiency, the b 
quark fragmentation, and the effects of the underlying event in defining the jet energy. The 
parameters in the smearing procedure are adjusted to account for these uncertainties, the 
input distribution is smeared, and the difference between the standard smeared distribution 
and the new smeared distribution is used to estimate the bin by bin systematic uncertainties. 
The uncertainties are added in quadrature to extract a total systematic uncertainty. Table || 
contains the bin by bin systematic uncertainties. 

6.2.4 b Jet p^ in Definition 

For future comparisons to theoretical predictions on overall normalization, we need to define 
a p™" - threshold for the recoiling b quark. The standard definition is to take the px value 
where >90% of all decays pass the kinematic cuts. In this case, we need to find the point 
where > 90% of all jets have Ex > 10 GeV. We begin with the b px spectrum shown in 
figure H and apply the resolution smearing to this distribution. We weight each bin in the 
Px spectrum by the probability that a b quark with that px would give a jet with Ex > 10 
GeV. Integrating the resulting weighted distribution gives a 90% p™ m value of 20.7 GeV/c 
for the b jet. 



28 



— p T Dis;ribution 




1 ~ 

_J I I I I LJ I I I I I I I I I I I I I I I I I Lj I I I I I I I lJ I I I I I I I I LlJ I I I L_ 

10 20 30 40 50 60 70 80 90 100 

p T (GeV/c) or E T (GeV) 



Figure 9: The best match p\ distribution, overlayed with the smeared distribution (dashed) 
and the data E T measurement. The process is reversed to take the data E T distribution to 
a pt distribution. 
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Figure 10: The differential E T cross section, for > 9 GeV/c, < 0.6, > 10 GeV, 
\r] b \ < 1.5, compared to theoretical predictions. The data points have a common systematic 
of ± 9.5%. The common uncertainty in the theory points comes from the muonic branching 
fraction and fragmentation model. The theory points do include uncertainties from the 
smearing procedure. 
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Variation 


10 — 15 GeV E T 


15 — 20 GeV E T 


20 — 25 GeV E T 


Energy Scale 
Underlying event 
Calorimeter Resolution 
Jet Finding 
b Fragmentation 


+ 7.2% - 4.6% 
+ 0.27o - 0.2/o 
+ 4.4% - 4.2% 

± 2.6% 

+ 1.0% 


+ 4.7% - 3.5% 

n 10/ n oO/ 
+ 0.17o - 0.27o 

+ 2.6% - 2.5% 

± 0.7% 

- 4.0% 


+ 9.1% - 7.3% 
+ 0.2/0 - 0.2/o 
+ 4.1% - 4.1% 

± 1.0% 

- 4.7% 


Total 


+ 8.9% - 6.7% 


+ 5.4% - 5.9% 


+ 10.0% - 9.6% 




25 — 30 GeV E T 


30 — 40 GeV E T 


40 — 50 GeV E T 


Energy Scale 
Underlying event 
Calorimeter Resolution 
Jet Finding 
b Fragmentation 


+ 12.5% -10.2% 
+ 0.3% - 0.37o 
+ 3.5% - 3.5% 

± 0.3% 

- 4.4% 


+ 16.5% -13.4% 
+ 0.4% - 0.4/o 
+ 0.9% - 0.9% 

± 0.0% 

- 3.4% 


+20.7% -16.5% 
+ 0.4/0 - 0.4/o 
+ 4.5% - 4.5% 

± 0.0% 

+ 1.6% 


Total 


+13.0% -11.6% 


+16.5% -13.9% 


+21.2% -17.2% 




50 — 100 GeV E T 






Energy Scale 
Underlying event 
Calorimeter Resolution 
Jet Finding 
b Fragmentation 


+27.8% -21.3% 
+ 0.4% - 0.4% 
+12.7% -12.7% 

± 0.0% 

+ 1.0% 






Total 


+30.6% -24.8% 







Table 8: Systematic uncertainties for each bin in the fi — b differential jet Et distribution. 
There are bin to bin correlations for each systematic variation. 

6.3 Comparison with NLO QCD 



In figure 10, we show a comparison of the differential jet Ex cross section, 



— > 9 GeV/c, \rf \ < 1, |r/ b | < 1.5, E^ > 10 GeV) 
dE^ 

to a prediction from the fi — b model discussed in section |3|. There is a 9.5% common uncer- 
tainty in the measured points, coming from the jet probability fit (5.8%), the b jet tracking 
efficiency (5%), the muon acceptance and identification efficiencies (3.9%), the luminosity 
normalization (3.6%), and the remaining background subtraction (2%). This common un- 
certainty is displayed separately. The uncertainty in the model prediction represents the 
uncertainty from the muonic branching fraction (5%) ||, the acceptance of the muon pp cut 
from variations in the fragmentation model (5%), which are common to all points, and the 
uncertainties associated with to Et smearing. The data has an integral value of 586. ± 



31 



61.8 pb, while the model predicts an integral value of 383.5 ± 5.9 pb. 
In figure [□], we show the unsmeared differential jet pt cross section, 



> 9 GeV/c, |^| < 0.6, \rf \ < 1.5) 

dp T 

compared to the b j?t prediction from the fi — b model, where we have included systematic 
uncertainties associated with the resolution smearing on the measured points. Again, the 
common normalization uncertainties are displayed separately. 

In figure O, we show a comparison of the differential — b) cross section, 



<]<J (p£ > 9 GeV/c, \rf \ < 0.6, E T > 10 GeV, \r] b \ < 1.5) 



to the predictions from the fi — b model. The uncertainty in the theoretical prediction 
represents the uncertainty in the muonic branching fraction and fragmentation model only. 

While we find qualitative agreement in shape between the measured distributions and 
model predictions, there are some differences. To investigate in more detail, we present in 
figure the experimental results minus the model prediction, scaled to the model prediction 
for the Et, Pt, and 5(f) distributions. The Et(pt) distributions have similar shapes for 
Et(pt) > 20 GeV( 35 GeV/c), but different normalizations. At lower values of Et (pt), 
the measurements and predictions are in agreement. The data 5(f) distribution is somewhat 
broader than the model predictions, with enhancement in the region tt/A to 37r/4, as well as 
being at consistently higher values. We have also shown how the model prediction changes 
with change of the renormalization and factorization scale, by plotting the prediction for scale 
/x /2 minus the prediction for /i , scaled to the prediction for /i . The integral cross section 
increases by 7%, with very little change as a function of Et or p^. In the 5(f) distribution, 
the /io/2 prediction is uniformly larger than the fio prediction, except for the region 5(f) ~ 7r. 

Recent work has shown that the addition of an intrinsic kick to a next-to-leading order 
QCD calculation improves the agreement between measurements and predictions for both 
direct photon production ]33J and charm production ||. We have investigated the effects of 



additional intrinsic fop i n the fx — b model. We use a Gaussian distribution with mean and 
adjustable width to model the magnitude of the kick, with a random azimuthal direction. 
With widths of 2 - 4 GeV/c, we find that the dominant effects occur for 5(f) < 1 radian. 
The cross section for 5(f) < 1 is predicted to change by approximately 7% with a width of 
4 GeV/c. With the current statistical uncertainties at small 5(f) (ranging from 25 - 100 %), 
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MNR Mode 




20 30 40 50 60 70 80 

Pt (GeV/ c) 

Figure 11: The differential p t cross section, for pj,>9 GeV/c, |ry M | < 0.6, |?7 b | < 1.5, 
compared to theoretical predictions. The data points have a common systematic of ± 9.5% 
and there are correlated systematic uncertainties. The uncertainty in the theory curves 
comes from the muonic branching fraction and fragmentation model. 
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Figure 12: The differential 5(f) cross sections, for pf^ > 9 GeV/c, |r/ M | < 0.6, > 10 GeV, 
\rf \ < 1.5 compared to theoretical predictions. The data points have a common systematic 
of ± 9.5%. The uncertainty in the theory curves comes from the muonic branching fraction 
and fragmentation model. 
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Figure 13: For the (a) da/dE T , (b) da/dpT, and (c) da/dScf) distributions, we plot the 
difference between the data measurement (filled circles) and the model prediction, scaled to 
the model prediction. There is a common systematic uncertainty of 9.5% in all the points, 
which has not been included in the error bar. The open circles are the model prediction 
for renormalization scale of /io/2 minus the model prediction for /i , scaled to the model 
prediction for /i . 
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we are unable to distinguish effects at that level. Similarly, the dominant effect in the p\ 
distribution occurs in regions where we have no sensitivity (p\ < 20 GeV/c). We conclude 
that the addition of intrinsic with width of 4 GeV/c does not account for the difference 
between the model prediction and the measurement. 



7 Conclusions 

We have presented results on the semi-differential \i — b cross sections as a function of the b jet 
transverse energy (da/dEx), b transverse momentum (da/dpr), and the azimuthal opening 
angle between the muon and the b jet (dcx/d£0). These results are based on precision track 
reconstruction in jets. The effects of detector response and resolution have been unfolded to 
translate the results from b jets to b quarks. We have compared these results to a model based 
on a full NLO QCD calculation |22| . We have investigated the effects an additional intrinsic 



&t and find that it cannot account for the difference between the measurements and the 
model prediction. Unlike previous CDF measurements |], |[], a normalization change alone 
does not account for the differences between this measurement and the model prediction. 
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